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^ • Abstract 

The major assumption of the Lorentz-Lorenz theory about uniformity of local fields and atomic 
polarization in dense material does not hold in finite groups of atoms, as we reported earlier [A. E. 
Q ; Kaplan and S. N. Volkov, Phys. Rev. Lett. 101, 133902 (2008)]. The uniformity is broken at 

Q^' sub-wavelength scale, where the system may exhibit strong stratification of local field and dipole 

o ■ 

. polarization, with the strata period being much shorter than the incident wavelength. In this paper, 
O , 

■ ^® further develop and advance that theory for the most fundamental case of one-dimensional 
>>■ 

r~| ' arrays, and study nanoscale excitation of so called "locsitons" and their standing waves (strata) 



X 



that result in size-related resonances and related large field enhancement in finite arrays of atoms. 
The locsitons may have a whole spectrum of spatial frequencies, ranging from long waves, to 



(N 
> 

XT) 

' an extent reminiscent of ferromagnetic domains, - to super-short waves, with neighboring atoms 

o ' 

' alternating their polarizations, which are reminiscent of antiferromagnetic spin patterns. Of great 

, interest is the new kind of "hybrid" modes of excitation, greatly departing from any magnetic 

I analogies. We also study differences between Ising-like near-neighbor approximation and the case 



where each atom interacts with all other atoms in the array. We find an infinite number of 
"exponential eigenmodes" in the lossless system in the latter case. At certain "magic" numbers 
of atoms in the array, the system may exhibit self-induced (but linear in the field) cancellation of 
resonant local- field suppression. We also studied nonlinear modes of locsitons and found optical 
bistability and hysteresis in an infinite array for the simplest modes. 

This paper is now published, with very minor changes, in Phys. Rev. A 79, 053834 (2009). 
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I. INTRODUCTION 

This paper is a theoretical extension of our recent Letter fll] on nanoscale stratification 
of local field and atomic dipole excitation in low- dimensional lattices driven by a laser at 
the frequency of the resonant atomic transition. We focus here on the most fundamental 
case of a one- dimensional (ID) array of resonant atoms, and construct a detailed theory 
of both linear and nonlinear interactions in the system, resulting in such phenomena as 
sub-wavelength spatial modulation (stratification) of polarization and local field, long-wave 
and short-wave stratification, size- related resonances and large field enhancement, "magic" 
numbers, ferromagnetic- and antiferromagneticlike atomic polarization, optical bistability 
and hysteresises, etc. In addition to the results l| and their derivations, we also present 
new results (i) on the size-related resonances using a many-body approximation involving 
interactions of each atom with all other atoms in the array, beyond the Ising-like approxi- 
mation whereby atoms only interact with their nearest neighbors, (ii) on traveling locsiton 
waves and their dissipation, as well as an estimate of the maximum size of the ID array 
to support most of the effects discussed here, (iii) on hysteresises and optical bistability in 
arbitrarily long arrays, and (iv) general mathematical consideration of dispersion relation in 
the self- interacting arrays for dipole approximation and beyond it. 

It is well known that optical properties of sufficiently dense materials are substantially 
affected by the near-field interactions between neighboring particles at the frequency of the 
incident field, in particular, quasi-static (non-radiative) dipole interactions. The best known 
manifestation of this fact is the local vs incident field phenomenon and related Lorentz- 
Lorenz or Clausius-Mossotti relations for dielectric constant as a nonlinear function of the 
number density. The microscopic (or local) field (LF), El, acting upon atoms or molecules 
becomes then different from both the applied and average macroscopic fields because of 
inter-particle interaction. In the most basic case, that relation between El and the average 
field E is El = (e + 2)E/3, where e is the dielectric constant of the material at the laser 
frequency u. It is worth noting that we are not interested here in the relation between e and 
the number density of the material, since in the case of small ID or 2D arrays the issue is 
moot. For the same reason, it also makes sense to deal directly with the incident field, Ejn, 
instead of the averaged field. In most theories of those interactions a traditional standard 
(and at that often implicit) assumption, refiected also in the above formula, is that the local 
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field and polarization are uniform in the near neighborhood of each particle at least at the 
distances shorter than the wavelength of light, A. This essentially amounts to the so called 
"mean-field approximation" . 

It was shown by us [1], however, that if the local uniformity is not presumed, then, 
under certain conditions on the particle density and their dipole strengths, the system of 
interacting particles is bound to exhibit periodic spatial variations of polarization and local- 
field amplitude. These variations result in sub-wavelength strata with a nanoscale period 
much shorter than A. Under certain conditions, the system may exhibit an ultimate non- 
uniformity, whereby each pair of neighboring atoms in ID arrays has their dipoles counter- 
oscillating with respect to each other, i. e. their excitations and thus local fields have opposite 
signs. 

To a certain extent this is reminiscent of the situation in magnetic materials with fer- 
romagnetic vs antiferromagnetic effects. Indeed, the mean-field approximation which is at 
the root of the Curie- Weiss theory of ferromagnetism jsl, is based on the assumption of 
uniform polarization of all the neighboring magnetic dipoles even without external magnetic 
field, when taking into account their interaction with each other. Contrary to that, the 
Ising theory [3], which does not make the assumption of uniformity, showed that even in 
the near-neighbor approximation for that interaction, the excitation may result in counter- 
polarization of neighboring atoms in ID arrays, and thus in a new phenomenon of antifer- 
romagnetism. Note here that in these effects there is no notion of "local" vs "external" field 
phenomenon: in a "pure" case of either of the magnetic effects, no external field is applied; 
the effects here are the result of self-organization of permanent, "hard", atomic dipoles with 
pre-existing dc dipole fields, without any "help from outside" . 

In this lies a profound difference between dc magnetic material phenomena (ferro- and 

n 

antiferromagnetism) on the one hand - and the effects considered by us here and in [Jj on 
the other hand, all of which are based on the optical (or, in general, any other quantum or 
classical resonance) excitation of atoms (or other small particles, e. g. quantum dots, clusters, 
small-particle plasmons, etc.). While magnetic dipoles in ferromagnetics are nonzero even in 
the absence of an external field (we may call them "hard" dipoles), the oscillating dipoles (in 
the linear case) can be induced only by the driving field at the near-resonant frequency, so 
they can be called "soft" dipoles; without such a driving their polarization vanishes. Because 
in dense material the atoms actually are acted upon by local field, the response of each one 
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of them may differ from the others by phase and amphtude (or even direction), with some 
of the dipoles fully suppressed while others fully excited. Thus the effects considered here 
are induced by the interplay of external and local fields, which put the entire phenomenon 
squarely into the domain of relations between the local and incident fields. Because of 
that, since the phenomenon depends strongly on the characteristics of the incident field 
(polarization, frequency, and, in the nonlinear case, intensity), the spatial modulation of 
the dipole excitation and local field can vary substantially. This results in a wealth of 
different patterns, some of them reminiscent of the ferromagnetic, other of antiferromagnetic, 
but the most of them forming all kinds of hybrid patterns. The complete cross-over from 
ferromagneticlike to antiferromagneticlike state of the system with all the intermediate states 
can be attained then by simply tuning laser frequency. 

Another significant difference here is that the system size is small. Provided there is 
sufficiently strong interparticle interaction, the new phenomenon can occur in the vicinity of 
boundaries, lattice defects, impurities or in sufficiently small group of atoms; recent advances 
in technology allow fabrication of nanoscale structures with small numbers of atoms. Thus, 
our theory emphasizes phenomena in relatively small ordered arrays of interacting atoms, 
in contrast to, e.g., microscopic models of ferromagnetism that mostly focus on averaged, 
"thermodynamic" perspective on sufficiently large systems. This brings forward a new set 
of nanoscale phenomena. Harking back to ferromagnetic systems, this new emphasis may 
reveal similar phenomena for nanoscale magnetic systems, which could be an exciting topic 
for a separate study. 

Our choice of ID and 2D dielectric systems based on the arrays or lattices of atoms, 
quantum dots, clusters, molecules, etc., allows to control anisotropy of near-field interaction. 
It also eliminates the issues of electromagnetic (EM) propagation being modified by the 
effects as the EM wave propagates through the structure (especially if it propagates normally 
to the lattice). 

If local uniformity is broken by any perturbation, the system may exhibit near-periodic 
spatial sub-A patterns (strata) of polarization. In general, two major modes of the strata 
transpire: short-wave (SW), with the period up to four interatomic spacings, la, and long- 
wave (LW) strata. The strata are standing waves of elementary LP excitations (called 
locsitons in pLj]) having a near- field, electrostatic, nature and low group velocity. 

In the first approximation, the phenomenon is linear in the driving field, and the locsitons 
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may be excited within a spectral band much broader than the atomic hnewidth. It can be 
viewed as a Rabi broadening of an atomic hne by interatomic interactions. The strata are 
controlled by laser polarization and the strength of atom coupling, Q, via atomic density, 
dipole moments, relaxation, and detuning. Once \Q\ > Qcr = 0{1), the LF uniformity can 
be broken by boundaries, impurities, vacancies in the lattice, etc. A striking manifestation 
of the effect is large field resonances due to locsiton eigenmodes in finite lattices, and - at 
certain, "magic" numbers of atoms in the lattice - almost complete cancellation of field sup- 
pression at the atomic resonance; saturation nonlinearity results in hysteresises and optical 
bistability. 

The paper is structured as follows. In Section [ITl we derive the main equations for self- 
interacting atomic lattices of arbitrary dimensions using two-level (nonlinear in general) 
model for atomic resonances and dipole-dipole interaction between atoms, while Section UTTl 
is on specific equations for linear infinite and finite ID arrays. In Section HVl we develop the 
general theory of locsitons and derive the dispersion relation. In Section |V] we study locsiton 
band formation, size-related resonances due to standing waves of locsitons (strata) and local- 
field enhancement. In Section lVll we concentrate on detailed theory of resonances beyond the 
near-neighbor approximation, including evanescent solutions (see also Section [IV]). Magic 
numbers are considered in Section IVIIi In Section IVIIIt we study the effects of losses 
on locsiton excitation, depth of penetration, and traveling locsiton waves. Section IIXI is 
on nonlinear locsiton modes, in particular, optical bistability and hysteresis. Section 1X1 
addresses potential applications of locsitons and their analogies in other physical systems. In 
Conclusions, Section IXll we summarize our results. Appendix |A] is on general mathematical 
aspects of dispersion relations for ID arrays. 

II. MAIN EQUATIONS 

Our model is based on the near-field dipole atomic interactions, with the incident fre- 
quency u being nearly resonant to an atomic transition with a dipole moment at the 
frequency cuq. In the linear case, i. e. when the laser intensity is significantly lower than that 
for the quantum transition saturation (see below), the result of this model coincides with 
that of a classical harmonic oscillator formed by an electron in a harmonic potential with 
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the same resonant frequency uq and with the dipole moment 

141 = ^^/X^, (2.1) 

where Ac = 2TTh/mc is the Compton wavelength of electron. 

In a standard LF situation, A ^ la, where A = 2ttc/uj, the field of an elementary dipole 
with the polarization p is dominated in its near vicinity by a non-radiative, quasi-static 



lis dominant term 
A Section 72). At 



(and only electric) component, which is strongly anisotropic in space, 
in the near-field area, |r' — r| = ro <^ A, attenuates as I/tq (see e.g. 
\da\ <^ X the amplitude of the field of an oscillating dipole located at r' induces a field 
at the point of observation r with the amplitude coinciding with that of an elementary static 
dipole with the same polarization p': 

„ , , , 3u(p' ■ u) — p' , , 

Edp(r',r) = ^i;,-^, (2.2) 

where u = (r — r') / is a unit vector in the direction of observation, and e is a background 
dielectric constant. 

We will model a resonant atomic transition by a basic two-level atom in a steady-state 
mode under the action of a field, E exp{— iut) /2 + c.c, with the amplitude E and frequency 
u, and assume that la ^ Ma|/e, so that the wave functions of neighboring atoms do not 
overlap. Using a semi-classical approach standard in LF theory of resonant atoms [s], ^, we 
can now find the atomic polarization as: 



2\da\^EAN 



(2.3) 



where p is the polarization amplitude [the full polarization is then pexp{—iujt)/2 + c.c.]; 
AA^ = Ni — N2 is the population difference, with Ni and N2 being atomic populations 
at respective ground and excited levels (^"1 + = 1), ^ = TAu = Tiu — ujq) is a 
dimensionless detuning from the resonant frequency ujq of the two- level atom, T = 2/r is a 
transverse relaxation time (the time of polarization relaxation) and F is the (homogeneous) 



linewidth of the linear resonance 



24| . In turn, the steady-state population difference is 



1 + 52 



AiV = iV^^ ( 1 + rr ' ) , (2.4) 



where r is a longitudinal relaxation time (life-time of the excited atom), and A^'''^ is an equi- 
librium population difference at the system temperature due to Boltzmann's distribution; 



in optics one can usually assume A^'^'^ ~ 1, so that 
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2/7712 n -1 



1 + 



1 + 52 



= fNd\E\ 



sat 



(2.5) 



where E^^^ is the saturation intensity, and /nl(I-EP) is nonlinearity due to saturation. Sub 
stituting fl2.5p into (12.31) . one obtains a closed- form expression for the polarization, 

214'' 



-E 



hT{5 + i 

For a classical harmonic (linear) oscillator we have: 

Ee2 



(2.6) 



P 



(2.7) 



uJoTm{6 + i) 

where m is the mass of electron. 

The local field EL(r) at each atom is the incident laser field, Ein, plus the sum of the 
near-fields, Edp(r',r) (12. 2p . induced by all the surrounding dipoles at r' acted upon by the 
respective local fields EL(r'), i.e. 



EL(r) = Ei,(r) + ^Edp(r',r) 

latt 



(2.8) 



latt 



where ^j^tt denotes summation over the entire lattice or array. To obtain a closed-form 
master equation, e. g. for EL(r) alone, we use (12.61) to write: 



EiYr) = E,,(r) 



Q 



E 

latt 



la 



X {3u[EL(r') ■ u] - EL(r')} /NL[|EL(r' 



'M21 



(2.9) 



where EL(r) are local fields only at the locations of atoms in the lattice and not at any other 
points inside or outside it; Q = Qa/{S + i) is a tuning-dependent strength of dipole-dipole 
interaction, and the maximum absolute strength, Qa, is 

8|4P 4:aXo{\da\/e)^uo 



Qa 



'■a 



(2.10) 
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where a = e^/hc ~ 1/137 is the fine-structure constant, and nonhnear factor /nl is as in 
(12. 5p . For a classical harmonic oscillator, (12.11) . we have 

where = e'^/rriec'^ ^ 2.8 x 10~^ nm is the classical radius of electron. 

Eqs. (12. 8p . (12.91) reflect many-body nature of the interaction. A conventional approach to 
local flelds within the Lorentz-Lorenz theory is to look for a self-consistent solution for the 
flelds in this interaction, with an assumption, however, that they are uniform (the mean- 
fleld theory), i.e. to set EL(r) = EL(r') and also use an encapsulating sphere around the 
observation point. These assumptions effectively shut out any strong spatial variations of 
the atomic excitations and local field that may exist at the inter-atomic scale. That is where 
we depart from the Lorentz-Lorenz theory; none of those assumptions are used here, and 
our approach is to use general expressions (12.81) or (12. 9p and seek straightforward solution 
for them. 

We will see below that the major critical condition for the phenomenon to exist and 
be observable at least at other optimal conditions is that Qa exceeds some critical value, 
Qa > 0(1). Three parameters are critical in this respect: atomic dipole moment da, the 
spacing between atoms, la, and the atomic linewidth, T, since Qa oc \da\'^/(Tll). To get an 
idea of whether the above critical condition is realistic, let us look first at the case of a gas-like 
collection of atoms, with the relaxed requirement on the spacing /„. Large dipole moments 
and narrow resonances in, e.g., alkali vapors or CO2 gas j^, in solids [7|, quantum wells 
and clusters may greatly enhance the phenomenon and allow for la from sub-nanometer to 
a few tens of nanometers. Considering an example with la ~ 100 A, corresponding to the 
volume density of ~ 10^^ cm~^, \da\/e ~ 1 A, Aq ~ 1 yum, e = 1, and T/uJo ~ 10^^, all 
of which are reasonable data, we obtain Qa ~ 10^, which provides a margin large enough 
to see all the effects discussed here. It is also of interest to roughly estimate what is the 
upper limit for Qa- To that end, consider the extreme situation of ~ \da\/e (solid-state- 
or liquid-like packing of participating atoms), in which case we have the ceiling for Qa as 

_ 8e=^ _4a Ap up 
"^^^^ ehVldal 7rei\da\/e)V ^''''^ 

Even taking into consideration significant line broadening, F, Qceii niay exceed unity by 
many order of magnitude, thus providing huge margin for the existence and observation of 
locsitons and related effects. 
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FIG. 1: ID array of atoms and laser light with different polarizations; light is incident normally to 
the plane of the graph. 

III. ID ARRAY OF ATOMS: LINEAR CASE 

We consider here the most basic model of a ID array of atoms lined up along the z 
axis, spaced by /„, and driven by a laser propagating normally to the array and having an 
arbitrary polarization (Fig. [I]). In the linear case, I-ElP -C -Egat; when in (12.91) /nl = 1, 
it is sufficient to consider effects caused by linearly polarized light with either one of the 
two mutually orthogonal polarizations; any other polarization (e. g. a circular one) can be 
treated as a linear combination of those two. In the case of a ID array, natural choices for 
these two basic configurations are: 

(a) the incident field Ein is parallel to the z axis, Ejn || (and the dipoles line up "head- 
to-tail"; we will call it || -configuration), and 

(b) Ein is normal to the z axis, Ein -L ("side-by-side" lineup; ±-configuration). 

The general solution will be a linear vectorial superposition of these two. This choice of 
the basic configurations is dictated by the simplicity of the resulting polarization of the 
local field. Indeed, in both cases, it follows that the polarization of local field is parallel 
to that of the incident field. El || Ein, so we can use scalar equations for all fields. Using 
the dimensionless notation Sn = [-EL(i"n)/-E'in](p), where (p) denotes polarization, (p) = || or 
(p) = _L, and recalling that now u || e^, we write ( 12. 9p for both configuration as 




1, if n = l,...,N, 



and = otherwise. 



(3.1) 



where is a form-factor due to polarization configuration: F\\ = 1 and F± = —1/2. 

If the ID array is infinite, —>■ oo, or sufficiently long, ^ 1, it is instructive to rewrite 
(13. ip in the form 

S^ + SQF,,, g |z^ = l' (3-2) 

where S = Yl'jLiJ^^ ~ 1.202057. The sums over \j — in (13.11) . (13. 2p converge rather 
fast, hence S* — 1 is not too large, see also Appendix |X] below. Eqs. (13.11) and (13. 2p . the same 
as master equation (12.91) . represent the case of fully-interacting arrays (FIA), whereby each 
atom "talks" to all the other atoms in the array, which presents a challenge to an analytical 
treatment. Of course, a linear equation (13. ip for Sn is solved analytically using a standard 
linear algebra approach with matrices. However, analyzing the results for large-size arrays, 
A^ ^ 1, in particular analytically finding all the resonances in (£^ri)max('^), can only be done 
by using numerical matrix solver, even if we neglect the dissipation. 

Thus, there is a need for a simple approximation that would preserve most of the qual- 
itative features of the phenomenon, yet could be easily analyzed analytically. This can be 
done by using the near-neighbor approximation (NNA), similar to that of the Ising model 
of (anti)ferromagnetism, in which the full sum in (13. ID or (13. 2p is replaced by the sum over 
the nearest neighbors, 

H ^^(^n-l + ^n+l) = 1; 

So = Sn+i = 0. (3.3) 

In the ultimate two-atom case, A^ = 2 |l|, the two approaches merge. The further two- 
near-neighbors approximation (2-NNA) and even three-near-neighbors approximation are 
considered in Section |VT] below. We found, however, that in general, a full summation (FIA) 
in (13. ip on the one hand, and NNA (13.30 as well as 2-NNA on the other hand, produce 
qualitatively similar results that differ by a factor of 0(1). 

Since effects discussed here are most pronounced in relatively small systems or in the 
small vicinity of perturbations in large lattices, it is natural to stipulate that the incident 
field within the array is uniform, unless otherwise is stated; however, this condition can 
readily be arranged even for an array larger than A. One of the solutions for the local field 
(and atomic excitation) in the infinite ID array (or a sufficiently long one, whereby we can 
neglect edge effects) is also uniform. We will call it the "Lorentz" solution, S, to be found 
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from (13. 2p by setting £n = £j = £ as: 

£ 1 - 1 

1 + Q5F(,) 1-[(5ll)(p)/(5 + ^)]' 

(5LL)(rt = -5g.F(,), (3.4) 

where we introduced polarization-dependent parameter which determines Lorentz- 
Lorenz shift at 6 = (see below); is a measure of the polarization-related strength 
of interaction. Eq. (13.41) may be viewed as a ID counterpart of the Lorentz-Lorenz relation 
for local field; notice, however, that the field £ is strongly anisotropic with respect to the 
polarization. The spectral behavior of \£\ is depicted in Fig. [2] with thicker dashed curves in 
all the graphs. If Qa 1, it shows, as one may expect, a deep dip at the atomic resonance 
frequency, i.e. at 5 = 0, 

l^(p)lmin = 1 I s:2 (3-5) 

(so that at the atomic resonance the local field is suppressed, as if it is pushed out of the 
array), and a strong new resonant peak appears at the shifted frequency 6 = 5ll due to the 
Lorentz-Lorenz effect: 



1 

\£{p) 



l^(p)lmax = 1 + (^LL = ic^ |2 > (3-6) 



mm 



whose nature is essentially similar to e. g. Lorentz-Lorenz resonance observed experimentally 

n 

in alkali vapors [6]. However, in the ID case considered here, the Lorentz-Lorenz shift and 
even its sign are polarization-dependent; in particular, 

(^Ll)iI =-Qa^, (5ll)± = (3.7) 

i. e. the Lorentz-Lorenz resonance is red-shifted for || -polarization of laser, and blue-shifted 
for ^-polarization. The Lorentz field £ and Lorentz-Lorenz shift in the near-neighbor ap- 
proximation, Eq. (13.31) . are determined by the same equations (I3.4l) - (l3.7p . where one has to 
set S = 1. 



IV. SPATIALLY-PERIODIC AND WAVE SOLUTIONS (LOCSITONS) 

We look now for solution of Eq. (13. 2p as the sum of uniform LP, £, Eq. (13. 4p . and 
oscillating ansatz A£ oc exp{±iqn), where q is an (unknown) wavenumber, similarly, e.g. to 
the phonon theory [8|, with the difference being that we have here an excitation of bound 



11 




to 



12 





1 


III/ 


100 




/■ 

1"^ /ft 


10 






1 












0.1 
0.01 


. (b) 

■ 


•ii' 

? 




100 200 



FIG. 2: Spectra of absolute maximum local field amplitudes near an atomic resonance {6 = 0) 
in the near-neighbor approximation with 5ll = 200: (a) locsiton resonances with 12 atoms in 
the array; (b) the same with 13 (magic number) atoms in the array; (c) merging and damping of 
locsiton resonances for a large number of atoms, N = 100. The curves show the amplitudes of local 
fields for size-related resonances, local fields due to Lorentz-Lorenz theory, and lower and upper 
amplitude envelopes of the resonances. 



electrons, and not atomic vibrations. Essentially, the locsitons may be classified as Frenkel 
excitons [8] because of their no-electron-exchange nature. 
The wavenumbers q are found via the dispersion relation 

cos(ng) 6 + i 



1 



n=l 



5i 



(4.1) 



LL 



12 




FIG. 3: Dispersion relations for an infinite array and negligible losses, i.e. normalized laser de- 
tuning the normalized wavenumber, g/vr. Curves show the dispersion for near-neighbor 
approximation, NNA (|4.2|) . two-near-neighbors approximation, 2-NNA (|6.6|) . and for fully inter- 
acting array (FIA), whereby each atom interacts with all the other atoms (j4.ip . which coincides 
with the analytical fit (jAlOj) . 

The behavior of D{q) in the lossless case, 5ll > ^ 1, is depicted in Fig. Oby the solid 
curve. Within NNA, we have to set S" = 1 and replace the sum in (14.11) by its first term: 

6 + i 



-DNNA(g) = cosg 



(4.2) 



->LL 



see Fig. [3l fine-dashed curve. Distinct oscillations emerge in the area between the two edges 
of the locsiton band. Their wavenumbers are determined from (14.11) or (14. 2 p by neglecting 
the dissipation, i. e. by assuming real q and 5^, ^ll ^ 1- band edges corresponds 

to the maximum (and positive) [-D(g)]max = 1 (same as for NNA, [-DNNA(Q')]max = -D(O) = 1), 
and thus to the Lorentz-Lorenz shift, 6 = The other, "anti-Lorentz" , edge is at the 
opposite side of the atomic resonance, and is determined by the minimum (and negative) 



Din). Thus we have: 



-'anti 



D(7l) 



which can be evaluated using an amazingly simple relation for the sum (14. ip . which is 
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actually valid for a more general sum and an arbitrary exponent p > 1: 

oo 

Ril[lPl^Il^ = _1 + J_ (44) 

n=l 

This can be readily proven by an appropriate rearrangement of the terms in the sums in (14.41) , 
see below flAllI) . In the case of p = 3 and infinite array, (14. ip . we have ^anti/i^LL = —3/4. 
Hence, the locsiton band is determined as: 

5 



5 3 , 
!> — >-- (or 
Oll 4 



2 



< 1 within NNA), (4.5) 



with well developed locsitons at 5ll ^ 1- Indeed, if the dissipation is neglected and the 
wavenumbers q are real, there is an infinite number of solutions for them within the limits 
(14. 5p : in this case the meaningful positive solutions, within the first Brillouin zone, are 
—71 < q < 71. The plots of the functions in the 1. h. p. of (14. ip and (14. 2 p are shown as 
vs q in Fig. [3l Based on (14.50 . the total width of the locsiton band in terms of 6 is thus 
(7/4)|(5ll|, if one accounts for the interactions of each atom with all the rest of atoms in the 
infinite array, whereas it is 2|5ll| in the near- neighbor approximation. 

To gauge the dipole-dipole interaction in the lattice, one can also introduce its Rabi- 
energy as: 

= -jr = < riuo. (4.6) 

..t brings about a locsiton energy band ~ 2hQR ^ tiT (if 5ll ^ 1) skin to those in solid-state 
8], photonic crystals {sJ], and electronic band-pass filters. 

One of the most interesting effects due to locsitons is a wide spectrum of the standing 
waves, strata, formed by them, see Fig. |H As we already noted, they can range from the 
very long ones, long-wave (LW) strata, with the maximum spatial period being double of 
the whole length of a ID array, similarly to the main mode of oscillation in a violin string, 
to short-wave (SW) strata whereby each dipole oscillates in counter-phase to each of its 
nearest neighbors. The latter mode is the strongest manifestation of the fact that the array 
of N atoms is a discrete-element resonator, akin a string of beads connected to each other, 
with the beads capable of the same kind of motion, when each individual bead is oscillating 
in counter-phase with its neighbor. Because in both cases the excitation is of dynamic 
nature and is due to external driving, the mechanical analogy of SW locsiton modes is more 
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FIG. 4: Strata patterns of excitation and local field in finite arrays, and their relations to the 
resonance tuning in the case of 64 atoms and 5ll = 200. Curves and patterns show long-wave, 
ferromagneticlike excitation near the Lorentz-Lorenz resonance (top curve in the left plot and top 
pattern in the center), counter-phase, antiferromagneticlike excitation near the anti-Lorentz edge 
of the band (middle curve and pattern), and hybrid excitation near the point of atomic resonance 
(bottom curve and pattern). Note that all the curves in the left plot show absolute normalized 
amplitudes of the local field. Since the fields are in general complex, their absolute amplitudes 
are positive, so the near-zero points in the schematic depiction of the "hybrid" mode actually 
correspond to the lowest points of the bottom curve in the left plot. 

adequate than that of a static ferromagnetic configuration vs antiferromagnetic one; see also 
Appendix |X] below. 

It can be immediately found, e. g. from consideration of the plots at Fig. HI upper curve, 
that the LW strata emerge at the laser tuning very near to the Lorentz-Lorenz resonance, 
i.e. in the limit 1 - < 1- In this case, in both FIA and NNA IK2h . the wave 

number ^lw and the respective spatial wavelength Alw 

27iL 



LW 



(4.7) 



In the infinite array, the longest A is up to 2'7r/a5LL, when 5ll ~ ^'^ ~ locsitons with 

longer wavelength get significantly suppressed. 
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The opposite limit, or the anti-Lorentz side of the locsiton band, 



3 S 

- + — < 1 (for NNA it is 1 + 



_5_ 

5ll 



«1), 



(4.8) 



defines SW locsitons, with 



~ 71, 



Asw/2 ~ /, 



(4.9) 



i.e. Asw/2 is the finest grain of locsiton structure, as one would expect from the "bead 
string" analogy. However, an incommensurability, i.e. mismatch between Asw/2 and the 
lattice spacing 1^, whose ratio is in general an irrational number, results in a strong spatial 
modulation of the SW, giving rise to a coarse LW-like structure. 

In the NNA, this coarse structure of SW mode has the half-period roughly the same 
as for a pure LW mode, ~ Alw This can be readily understood in terms of "beating" 
between the locsiton wavelength, 2iT/q, and the spatial scale of the discrete structure of 
the system, which is the normalized spacing between atoms, 1. Indeed, near the anti- 
Lorentz edge, 6 + 5ll ^ ^ll, we can write for a SW wavenumber: gsw = vr — Aq, with 
|Ag| ^ a/1 — (5/5ll)^ = Q'lW) and find the spatial oscillations as 



with Aq ^ gLW) which shows alternating, counter-phase motion of the neighboring atoms, 
(— 1)", modulated by a slow envelope, cos(?T,gLw)- Both the fine grain and the coarse mod- 
ulation may be well pronounced (Fig. HJ middle curve). At g = vr, the LW and SW periods 
converge to the same scale, Ala, see ( 14.12^ below. Using the phonon analogy, the LW loc- 
sitons may be viewed to a certain extent as counterparts to acoustic, and SW - to optical 
phonons. 

Between those limiting points - LW locsitons at the Lorentz end of the locsiton spectrum, 
S ~ ^LL, and SW locsitons at the anti-Lorentz end, 6 ~ — ^ there are all kinds of loc- 
sitons making chaotic looking strata (due to the above mentioned incommensurability, i. e. 
irrational ratio between a locsiton wavelength. A, and the array spacing, /„). However, same 
as in the chaotic motion, there are small islands of well ordered wave patterns, located at the 
spectral points where the ratio A/ la (or q/if) is a rational number, provided that the system 
is at a resonance, see (15. 2 p and (15.31) below. One can think of them as sort of hybrids of fer- 
romagneticlike and antiferromagneticlike behavior. Indeed, the purely antiferromagneticlike 



Sn oc cos(ngsw) = cos(n7r — nAq) 



(-l)"cos(nAg) 



(4.10) 
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SW locsiton at g = TT is formed by the atoms with alternating polarizations, 



■ ■ ■ T i T i T i • ■ ■ for ± polarization, and 

•■■ ^< >< >^ ■■■ for II polarization. (4-11) 

This happens at 5 ~ — 5ll within NNA and 5 ~ — (3/4)5ll for FIA. The simplest hybrid 
pattern is formed at g = 7r/2 [which is at 5 = within NNA, and = —3/32 for a fully 

interacting array (14.11) ]. whereby each second atom is non-excited, while the other atoms 
alternate their polarization: 

••• folofolofol--- for± polarization, and 
■ •■^o^o^o^o^o^--- (4.12) 

for II polarization. 

Examples of other simplest hybrid states within NNA are: q = 7r/3, q = 7r/4, etc.; in 
general, periodic patterns exist if g/vr is a rational number, i.e. for all the eigen-resonances 
in finite ID arrays within NNA, see Eq. (15.31) below. 

Finally, it is worth noting that the strata, albeit fast-decaying, exist even beyond the 
locsiton band, (14.51) . They are not, however, propagating locsitons, and their amplitudes 
exponentially decay with the distance. Consider the simplest case, NNA, with the losses 
negligibly small, 5^ > 5ll ^ 1 in (14.21) . In this case, cos^ g > 1 in (14. 2p . which indicates that 
the wavenumber g must be complex. Indeed, writing g = —ix for the "Lorentz" end of the 
band, > 1, and g = —ix + vr for the anti-Lorentz end, 6/6i, < —1 (see also Section IVll 
below), we have the solution for x as 



X = In [\5/5ll\ ± ViS/SLL^ - ij (4.13) 

and for the local fields as 

(^„)LLCxe^" and (^„)antiLL « e>^"(-l)" (4.14) 

for the Lorentz and anti-Lorentz ends of the band, respectively. Here, in the case of a 
semi-infinite array, the sign in (I4.13P has to be chosen such that the amplitude of the field 
vanishes at \n\ oo. These modes can be viewed as evanescent waves/locsitons, see also 
below, Section |VT1 One can note though that they still bear the signature of the respective 
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locsitons on either side of the band: long-wave, almost synchronous oscillations on the 
Lorentz end, and short-wave, phase-alternating oscillations on the anti-Lorentz end. In 
general, the same patterns hold in the FIA case. 

V. QUANTIZATION OF LOCSITONS AND FIELD RESONANCES IN FINITE 
ID ARRAYS 

Due to boundary conditions in (13.11) [or (13.31) within NNA], the array of atoms is a 
discrete-element resonator. This should result in locsiton quantization within the locsiton 
band (14. 5 p and corresponding size-related resonances of the local field. In a ID array with 
atoms, we have N coupled oscillators with the same individual atomic resonance frequency, 
and therefore, one should expect the original atomic line to be split into lines at most, 
with the collective broadened band (14.51) being replaced with those lines. Of course, when 
the dissipation (or finite linewidth of each individual line) is taken into account, those split 
lines will merge into one continuous locsiton band (14. 5 p if the array is sufficiently large, 

> ^I^llI (or > 2|5ll| within NNA). (5.1) 

The simplest result for the resonant line positions is obtained within NNA. Using the 
boundary condition (13.30 . i.e. Sq = £n+i = 0, we find that the longest locsiton half-wave, 
corresponding to the fundamental mode, is 

Ai/2 = (AT+l)/,, q^ = Tx/{N + l). (5.2) 

Thus, Ai/2 is the distance between nodes where LF zeros out, whereas the wavelengths A^ 
of eigen- locsitons and their eigen-frequencies 5^, with the quantum number 1 < k < N , are 
respectively 

Afc = Ai//c, = (^LLCOsgfc, 

qu = nk/{N + l). (5.3) 

[Note that the first Eq. (7) in fl)], which corresponds to the second Eq. (15.30 here, contained 
a typo (an extraneous vr in the cos argument) which we corrected here.] From these, only 
the resonances with odd k will be realized for a symmetric driving laser profile, in particular, 
the uniform one, E^^ = const (which is the most common case here), and with even k - for 
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an anti-symmetric one, say, E'in = const ■ (A^ + 1 — 2n)/{N — 1), where n is the sequence 
number of an atom in the array. In all other full set of resonances will be realized. 

In essence, the size-related locsiton resonances in discrete arrays are, to a limited extent, 
similar to any eigen-resonances in regular continuous, i. e. non-discrete, ID system. Examples 
can be found both in classical setting, e. g. a violin string, a Fabri-Perot resonator (as e. g. in 
a laser), and in quantum mechanics, from the resonances in a quantum well with infinitely- 
high walls, to electron gas in a fin.e .aye. Q, eleCo. in long .nolec.es fl, etc. TKe 
major difference here is that the number of eigenmodes, or resonances, in an array with 
elements is limited to N, in contrast to the theoretically infinite number of eigenmodes in 
continuous finite ID systems. 

The resonances for uniform driving within NNA are shown in Fig. [2] for = 200 in the 
case of AT = 12 [Fig. [2](a)], A^ = 13 [Fig. [2](b)], and A^ = 100 [Fig. EJ^c)]. One can readily 
find out that the lower amplitude envelope is 

S,,,U6) ^ 2\Sl (5.4) 
while the upper envelope of the resonant peaks within NNA for a uniform driving is 

\^\S\{ns + nJ^), if < 1, 
I = 2 1^1 otherwise, 

where n5 = (A^+ 1)/(2 v^^ll " '^^)- As A^ increases, the resonances merge and are suppressed 
at A^ = |5ll|0(1), see e.g. Fig. [2](c). However, even then £^i^ still exceeds the uniform field 
1^1 (13. 4p by a factor of 2. For A^ = 3A; — 1 (fc is a natural number), LF amplitude dips below 
the lower envelope |£^max|iow at 5 = — 5ll/2. At that frequency, within NNA, cosg' = —0.5, 
q' = 2n/3, and the SW period A = 3/^ is an integer of the atomic spacing, so only fine 
SW structure remains, resulting in an anti-resonance and in the strongest inhibition of the 
locsiton. 



VI. ID ARRAYS BEYOND THE NEAR-NEIGHBOR APPROXIMATION 

As we mentioned above, while the quantization of locsitons in finite arrays can be readily 
analyzed analytically within NNA, see the previous Sections IIIII and |Vl the situation with 
fully-interacting arrays (FIA) presents a challenge for an analytical treatment. 

Let us briefly outline general analytical and numerical results obtained so far. 



19 



A FIA locsiton band is not symmetric with respect to the atomic resonance, 5 = 0; 
it is shorter by the factor 3/4 on the anti-Lorentz side, see (14. 5p . as opposite to the 
NNA. 



• Respectively, FIA resonances are grouped tighter on the anti-Lorentz side of the band, 
albeit their number is the same as for NNA. Near the Lorentz side of the band, the 
NNA-predicted resonances coincide more closely with those obtained by FIA numerical 
calculations. 

The major source of these effects is the first factor, i. e. a strongly asymmetric (with 
respect to the detuning frequency, 6) shape of the dispersion relation. Fig. [31 A more 
detailed mathematical consideration of this problem, including a very good analytical fit for 
the dispersion relation is found in Appendix 1X1 below. Less significant, although interesting 
as far as the eigenmodes of ID arrays are concerned, is the fact that simple NNA eigen- 
wavenumbers (15. 3p . qk = 7Tk/{N + 1), obtained based on the NNA boundary conditions 
So = Sn+1 = [Eq. (13. Sp ] are not exact anymore. One has to use now more extended, 
"beyond-the-boundary" conditions (13.11) . which are the signature of FIA, whereby 



instead of just two end-points in NNA. To explore the problem, let us simplify it first 
by considering only non-dissipating atoms, > ^ 1, and rewrite the full- interaction 
dispersion relation (14. ip for this case asi 



with the boundary conditions (16.11) . 

The plot D{q) for real g's due to the dispersion relation (16. 2p is shown in Fig. [31 with the 
solid curve. The new qualitative difference now between FIA, (16.21) . and the NNA case, (14. 2 p 
without dissipation, i.e. cosg = D, is as follows. With D'^ < 1, the NNA equation has only 
real solutions for q, whereas Eq. (16. 2p . even within the locsiton band, — 3/4 < < 1, aside 



Sn = for all n < 1 and n > N, 



(6.1) 




(6.2) 



assuming D{q) real. The equation for the field is written then as: 




(6.3) 
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from real solution for q = q, has, as one can show, an infinite number of complex solutions, 
q, for each single real solution q (i.e. for the same D). All of them, in addition to fast 
spatial oscillation terms, have a rapidly rising/ falling exponential factor. These exponential 
modes are negligibly small almost over the entire array length, if ^ 1, and they need to 
be accounted for only very near the end-points of the array, where they are instrumental in 
zeroing out the field and excitation at the points n < 1 and n > N. 

Let us illustrate the formation of those exponential (or evanescent) modes and their role 
in boundary conditions for the iwo-near-neighbors approximation, 2-NNA, whereby the field 
Eq. (16.31) becomes: 



9 

£ - — 
" AD 



{£n-l + £n+l) + g(^n-2 + ^n+2) 



(6.4) 



with the boundary conditions for two pairs of end-points: 



£n = at n = 0,-l and n = N + 1, N + 2. 



(6.5) 



The dispersion relations approximating ( \6.2\\ will read now as: 

cos(2g) 
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cosg -|- 



Diq), 



see Fig. [3l long-dashed curve, and the locsiton band is determined by 



(6.6) 



7/9 < D{q) < 1. 



(6.7) 



The real solutions q of ( 16.60 are those for which cos^ q ^ they give rise to strata modes of 
Eq. (EH), 

£n oc exp{±iqn). (6.8) 

However, having in mind that cos(2g) = 2 cos^ q ^ ^, one can readily see that (16.61) has also 
solutions with cos^g > 1, i.e. those that correspond to exponential modes, with complex 
q = '?evn2- Introducing for those modes 



gevn2 = -iX + T^ 



(6.9) 



with real we obtain from (16.61) that for each given real q the exponent x is determined 
by: 

cosh X = 4 + cos g, (6.10) 
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and the respective exponential mode is 

(^„)2 = e±(^+'")" = e±^"(-l)". (6.11) 

Thus, these 2-NNA modes are antiferromagnetichke strata, modulated by fast exponents. 
Indeed, since 3 < coshx < 5, we have 1.76 < x < 2.3. 

Since the exponential, or evanescent, modes have so short "tails", they produce relatively 
small correction for the respective eigen- wavelengths, A„, and for g^, compared to the NNA 
oscillatory modes, (15. 2p . (15. 3p . so that we can look for the A„ corrected for 2-NNA at the 
points of resonances as: 

Afc N +1 + Ak , , 

» = (6.12) 

with the correction > 0, similarly to, e.g., oscillations in a violin string with a "soft" 
suspension at its ends. To find of the k-th resonance, we seek a solution for Sn as a sum 
of two modes: oscillatory and exponential ones. Then, for symmetric modes, i. e. with odd 
k, we have a full solution written as 



(i^n)odd = cos(gfcn) + Ccosh[x(gfc)?^](-l)", 

with n = n ^— , (6.13) 

where C is a constant. For anti-symmetric modes, with even k, we have: 

(^n)even = sin(gfcn) + Csinh[x(gfe)n](-1)". (6.14) 

Using now conditions (16.51) . we can write for the points n = —1 and n = 0, respectively, the 
following equations for the symmetric modes (I6.13p : 

Vfc A^ + 3 



(^^-l)odd = COS 



2 A^ + l + A 



Ccosh ( x^^^ ) = 0, (6.15) 



(^^o)odd = cos 



2 

-Kk N + 1 

TivTTTA 



+ Ccosh ( x^^^ ) = 0. (6.16) 
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For the anti-symmetric modes fl6.14p . one has to replace cos and cosh functions in (16.151) . 
(I6.16P with sin and sinh functions, respectively. From these two equations we can compute 
Afc and C. Indeed, approximating cosh(^) ^ sinh(^) in (I6.15P and (I6.16P or in the 

respective equations for anti-symmetric modes, since x(A^+ 1) ^ 1, we obtain two equations 
for Afc and C, which are readily solved for A^, and C for both symmetric and anti-symmetric 
modes as: 



2 -1 
■ tan 



sin(g 



NNAJ 



(6.17) 



gNNA [e^ + cos(gNNA) 

with gNNA = 

and C = gNNAA,(-l)^/2e->^(^+i), (6.18) 

where from (I6.10p . x ~ cosh^^[4 + cos(gNNA)], while A; = A; + 1 for symmetric modes, and 
k = k for the anti-symmetric ones. It is worth noting that the boundary conditions (16.51) at 
n = A^ + l,A^ + 2 are now satisfied automatically because of our choice of the coordinate n 
in (I6.13P and (I6.14p . For LW resonances, k <ti N, (16.170 reduces to 

Afc ^ ^ A,„ax ~ ^, (6.19) 

which is the maximum magnitude of A^, whereas for SW resonances, A^ + 1 — A;^A^ + 1, 
(I6.17P reduces to 

which is substantially smaller than the LW correction (I6.19p . One can readily see that A^ is 
a monotonically decreasing function of k. In the middle of the locsiton band, k ~ (A^+ l)/2, 
we have 

A, ^ ^ ^ ^. (6.21) 

vre^ In 

Interestingly, a fairly good fit to (I6.17P is provided by a much simpler formula: 



Afc ~ Ajjiax 



1 



k 



N + 1 



(6.22) 



Now, once the correction A^ is found, one can substitute q = % = 7rA;/(A^ + 1 + A^), (16.121) . 
into the dispersion relation (16. 6p and calculate the respective frequency detuning for the fc-th 
resonance, = Sk/Si^, with k = 1 being the closest to the Lorentz-Lorenz resonance, i.e. 
a LW mode, and k = N closest to the anti-Lorentz edge of the locsiton band, a SW mode. 
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A way to extend the 2-NNA approximation to a full-array interaction is to apply the 
result fl6.17p for the correction of the eigen-wavenumber, g^, but use it now in the full-blown 
dispersion relation fl6.2p . instead of the 2-NNA relation (16.61) . to calculate the resonance 
detuning 5k- The other avenue, of course, is to seek for higher order approximations. For 
example, one can take into account two edge sets of 3 points each, i.e. similarly to (16.51) . 
stipulate 3-NNA: 

8n = Q at n = 0, -1, -2 

and n = + 1, A^ + 2, + (6.23) 

Following the same route as for 2-NNA case, we obtain now, similarly to (16.101) . a second 
order equation for the exponential eigenmode complex wavenumber, q = gevns, for any given 
oscillation wavenumber, q, as 

16 cos^ q + cos q (16 cos q + 27) 

+ (104 + 16cos^g + 27cosg) = (6.24) 

[we leave it untransformed to coshx, unlike (16.91) . (16.101) . since in 3-NNA case the solution 
becomes more complicated than (16.111) . and that transformation does not simplify the prob- 
lem.] The higher-order equation for the exponential eigenmodes ensue the choice of a higher 
number of end-points in the boundary conditions. 

It has to be noted, however, that for specifying parameters of oscillating eigenmodes, the 
increase of the precision by accounting for higher-order exponential modes is of very limited, 
if purely academic, significance. Essentially, for large arrays, ^ 1, the combination of 
the NNA for predicting the wavenumbers q, (16.121) or (15. 3p . with these g's used then in 
(16. 2p to calculate the resonance eigen-frequencies 6, does already a good job. A further step 
in increasing the precision provided by the 2-NNA corrections to eigen-wavenumbers, is to 
again use the corrected g's in (16.20 for the same purpose, and it is more than sufficient. A 
special case is a small array, e. g. = 2,3,4, whereby it is actually preferable to simply 
solve the general ID-array equations analytically, similarly to [I] in the case of A^ = 2. 

VII. MAGIC NUMBERS 

A fundamental effect of self-induced cancellation of local-field suppression emerges near 
the atomic resonance, 5 = 0, at certain "magic" numbers A^. If ^ 1, the uniform 
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(Lorentz-Lorenz) LF (13.40 at 5 = is very small, \£\ ~ ^ [Eq. (13.50 ]. However, in the 

near-neighbor approximation, if 

A^=/cmmag + l, m = 1,2,3,..., (7.1) 

where m^ag = 4 is a "magic" number within NNA, locsitons at (5 = show canceled LF 
suppression at some atoms. The highest cancellation is attained at 5 = and = 5, 
with the atomic dipoles lining up as in (14.121) . whereby the LF amplitude at odd- numbered 
atoms is at maximum, |£^max| ~ 1/3, and the enhancement |i^max/^lonh ~ could be a 

few orders of magnitude; the LF at the two other atoms almost zeros out. 

The self-induced cancellation effect is produced by a standing wave with the nodes at 
atoms with even numbers. This results in a "virtual" size-related resonance at 5 ^ (i. e. at 
the exact atomic resonance), which manifests itself in the enhancement (the resonant peak 
transpires in |£^max/^| vs 5). Thus, the nature of magic numbers is the coincidence of the 
atomic resonance with one of the size-related locsiton resonances. 

The effect holds also for the interaction of each atom with all other atoms (13. ip : the 
magic number in (17. ip takes on a "devilish" likeness here: m^^g = 13. It is due to the fact 
that for 6 = and 5ll ^ 1, the first root q' of Eq. (14.10 with zero r. h. s., 

cos(ng') 



5:^^ = 0. (7.2) 



n=l 



has its property of g'/vr almost coinciding with a rational number, q'/n ^ 6/13 (13g'/67r = 
1.00026 . . .), see Appendix 1X1 below. The locsiton wavelength is A = 2iT/q' = (13/3)/^, and 
the lowest integer of A/2 to exactly match an integer of la is 13/^; which requires 14 atoms. 
We have now |£^max| ~ 2/15, with enhancement ~ 45ll/225. 

As we have shown in [1], 2D lattices can also exhibit ^^magic shapes" with similar proper- 
ties; the simplest one within NNA is a six-point star with an atom at its center, thus making 
the total number of atoms again A^ = 13. More details on magic shapes for 2D lattices will 
be discussed by us elsewhere. 



VIII. TRAVELING LOCSITON WAVES: VELOCITY AND PENETRATION 
DEPTH 

If the locsitons are waves, can they be excited outside the driving field area? and thus 
travel away from the Lorentz-Lorenz uniform local field area? how far can they travel before 
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being extinguished? how fast do they propagate? having in mind their dissipation, what is 
the size of a finite array to have well pronounced standing waves and size-related resonances? 

Of course, the locsitons can propagate in an array even in the areas unaccessible for the 
external (laser) field. If the spatial profile gradient of the driving wave is large enough, 
a LF excitation can be found beyond the driving field area; the terminological irony here 
is that the /oca/-field phenomenon is due to a nonlocal interaction, and the locsitons can 
propagate away from their origination point. This may happen when the external laser field 
is non-uniform and has a large gradient, e.g. when one entirely screens out a part of the 
array by imposing a "sharp knife" over the array. 

Even simpler and more transparent case is when only an end-point of a semi-infinite ID 
array [say, with n = 1 in (13.11) ] is illuminated by a laser field via a pin-hole. Eq. (13. ip has then 
non-zero (unity) r. h. p. for n = 1 only, and the r. h. p. zeroes out at all the other points. In 
this case, no Lorentz-Lorenz local field exists for any atom at n > 1, and the only field and 
atomic excitation passed along the array of atoms will be locsitons. This may be the best 
way to excite and observe "pure" locsitons, with them not being masked by any external, 
averaged, mean, etc., fields. With such an arrangement, the ID array (or a sufficiently thin 
atomic "cylinder" or carbon nanotube) becomes a true and effective waveguide for locsitons, 
capable of transmitting non-diffracting radiation and atomic excitation from one location 
(e. g. in opto-electronic circuits) to another. 

The main issue here is how far the locsiton can propagate. One can investigate it by study- 
ing the dispersion relation (14.11) or (14. 2p . which predicts not only the locsiton wavenumbers, 
q' = Re(g), for any given frequency detuning 6, but also their dissipation depth or distance 
(in terms of numbers of atoms) for each wavenumber, as N^is = 1/q" where q" = Im(g). 
Using only the NNA dispersion relation (14. 2p . amply sufficient here, since the calculation of 
the dissipation distance doesn't require the same precision as for the real wavenumbers, we 
find out that the exact solutions for q' and q" are determined by: 
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and 

1 



smh\q") = -- ( 1 



2V* 61^ 



'1 + 1 



From (18. ip . for very small dissipation, i. e. when 5ll > ^ 1, we have, as expected, 

cos(g') ^ — . (8.3) 

The dissipation length in terms of the numbers of atoms, A^^dis = ^/q", is readily found from 
(18. 2p . At the exact atomic resonance, 5 = 0, if 5ll ^ 1; have: 

<l"-T^, iVdis-I^LLl- (8.4) 

I^llI 

In general, in the most part of the locsiton band (except for the edge areas, 5ll ~ ~ 1); 
we find that 

A^dis ^ ^feh^'. (8.5) 



q" 



Finally, at the locsiton band edges, defined as — 5^ = 1, (18.21) results in 

q" ~ iVdi3 ^ (8.6) 

V^LL 

and it remains roughly the same as 5^ reaches 

Thus, in the most of the locsiton band, the dissipation length, in terms of the numbers 
of atoms, is A^dis = 0(5ll); which also determines the maximum size of array to still enable 
the size-related resonances, in agreement with Section |Vl 

Let us address now the characteristic velocities of the locsitons. Neglecting decay in (14. 2 p 
(5ll > ^ 1)) the group velocity of locsitons, Vgr = la{,duj/dq), is found as: 

la dS S 
Vgr = --!-, cosg = — , 8.7 
Tdq 6ll 



hence 



dq _ 1 

ds - "VIT^^ 



laJnl-Auj^ = vnJ6l^-6\ (8.8) 
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where Qr is the Rabi frequency of the self-interacting array (14.61) . and vr = fl-Rla is a 
characteristic Rabi speed, 

Oil) < 1, (8.9) 



— = a 

c 



ela 

which could be even slower than a typical speed of sound in condensed matter. This effect 
can be used e. g. for developing nano-size delay lines in molecular computers and in optical 
gyroscopes. The LW-locsiton phase velocity is 

Vpt ~ — . (8.10) 

IX. NONLINEAR EXCITATION OF A ID ARRAY; OPTICAL BISTABILITY 
AND HYSTERESIS 

So far we studied linear (in field) excitations of atomic ID arrays. The nonlinear in- 
teractions open door to a huge new landscape of effects. The simplest and very generic 
nonlinearity in a two-level system is the saturation of its absorption (12.41) . (12.51) . which 
translates into a nonlinear response of the atom polarization to the local field (12. 6p : this, in 
turn, nonlinearly affects the strength of interaction between atoms (12.91) . This represents a 
rare case whereby the nonlinear change (decrease) of absorption directly affects the eigen- 
frequencies of the system (ID array), by directly reducing the interaction. Many nonlinear 
effects are brought up in short pulse modes, e. g. discrete solitons, to be considered by us 
elsewhere. However, spectacular effects, such as hysteresis and optical bistability, emerge 
even in cw mode. 

To write a set of nonlinear equations for an infinite ID array, we first scale all fields to 
the characteristic saturation field, i^sat (12. 5p . instead of scaling to the incident field, E[^, so 
that the dimensionless local fields at the n-th atom, Yn, and the incident field, X, are 

X = Ei^/Esa.t, yn = En/Esa.t, (9.1) 

and the nonlinear counterpart of (13. ip and (13. 2p for the array is written now as: 

n-^..(^-.)E ,^.^„,3(^f,^^|,.|,, ^A-. (9.2) 



latt 



where is defined by (13.40 . so it covers either incident polarization. The nonlinear coun- 
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terpart of the NNA equation (13.30 is now 



Yn- 



2 vi + 52 + i>; 



n-1 



12 



None of this is an easy object even for numerical solution, let alone analytical one. We 
have, however, derived in [ij] a closed analytical solution for the nonlinear mode in the most 
fundamental system - an array of just two atoms - and found bistability and hysteresises in 
such a mode. In this paper, more as a matter of illustration, we find an analytical multistable 
solution and hysteresis for the simplest case of Lorentz-Lorenz uniform field; however, our 
computer simulations showed that multistability and hysteresises exist in the vicinity of each 
size-related resonance in the system. 

For the sake of simplicity, we consider here the case of Lorentz-Lorenz uniform mode. 
In the near-neighbor approximation (19.31) [the FIA results will essentially differ only by a 
factor of 0(1)], we have Yn = Yn-i = = Y, and thus the nonlinear equation for the 

uniform local field, Y, is 



Y 

or for the field intensity, |yp = y, 



1 - 



1 + 52 + |F|2 



X, (9.4) 



2 



It can be readily seen that the strongest nonlinear effect emerges near the Lorentz-Lorenz 
resonance, 6 ~ Assuming small losses, 5, 5ll ^ 1, we can analyze the threshold of the 
multistabihty and hysteresis mode by stipulating that y <^ ^ = '^ll — 5 ^ and 
thus (19. 5p can be further simplified as 

y[iy~5^^Ar + 6l^]/6t^^X'. (9.6) 

The threshold for the multistable solution y{X) of this equation is determined by the con- 
dition dX/dy = (ff'Xjdy^ = 0, which results in the critical requirement that 

A = 5ll - 5 > Ae, = VS, 

y = > y,, = A^LL, (9.7) 
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FIG. 5: Multistability of Lorentz-Lorenz mode in a long array: local Lorentz field amplitude Y vs 
detuning 6 for different normalized driving amplitudes X = £'in/i?sat for the case of 5ll = 100. 

Amazingly, as one recalls that X = Ei^/ Es^t and 5ll ^ 1, the critical (threshold) driving 
intensity E^^ to initiate multistability and hysteresis could be orders of magnitude lower than 
the saturation one, E^^^, which is mostly due to resonant nature of the effect that emerge 
in the vicinity of the Lorentz-Lorenz resonance. Thus, the required saturation nonlinearity 
is indeed just a slight perturbation to the resonant linear mode. This should be of no big 
surprise, since the nature of the effect is the same as in many other so called vibrational 



hysteresises 



i 



in resonant nonlinear systems, from pendulum [12l] to electronic circuits 



131], to optical bistability in a Fabri-Perot resonator IJ], to a cyclotron resonance of a 
slightly relativistic electron |15i]. In all these nonlinear resonances, it is enough for a narrow 
resonant curve to be tilted beyond its resonance width to reach multistability by becoming 
a "Pisa-like tower". 

The formation of bistability and hysteresis is depicted in Fig. [5], where the resonant curve 
Y{6) is shown for different driving X for the full equation (19.51) . The formation of a tri- 
stable solution results in the so called bistability, whereby only the states with the maximum 
and minimum intensities are stable, whereas the middle branch of the solution is absolutely 
unstable. To an extent, it is similar to one of the hysteretic patterns for the local field found 

n it 

for the case of the local field pUj and scattered light pj| of just two atoms; the hysteretic 
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resonance here is produced by a Lorentz-like, or ferromagneticlike, excitation of atoms with 
LW locsitons. This effect is also reminiscent of optical intrinsic bistability for uniform LF 
in dense materials predicted earlier in (l7] and observed experimentally in [l^. 

The above result for hysteresis and bistability was found for the uniform, Lorentz, field, 
i. e. near the Lorentz-Lorenz resonance, using a simple field distribution. However, it is 
clear from the above argument with a "Pisa-like tower", a similar hysteretic effect can be 
expected in the close vicinity of every size-related locsiton resonance within the locsiton 
band. Some of the most interesting hysteresises, though, are excited for the SW locsitons, 
when the nearest atoms counter-oscillate in an antiferromagneticlike fashion, as shown in 
|l] for the anti-Lorentz band edge for two atoms. This hysteresis results in a "split-fork" 
bistability and will be discussed by us in application to long arrays elsewhere. 



X. ANALOGIES TO LOCSITONS IN OTHER PHYSICAL SYSTEMS AND GEN- 
ERAL DISCUSSION 

It is only natural to expect the effects studied here to manifest themselves in other finite 
systems of discrete resonant oscillators that are externally driven and interact with each 
other. On the one hand, this may help to demonstrate and verify some major results reported 
here in much simpler and easily handleable settings; on the other hand, this may bring up 
new features in these other systems that escaped attention in seemingly well developed and 
researched fields. We consider here a few such systems. 

Perhaps, the most classical example would be a mechanical finite array of identical phys- 
ical pendulums with the same individual resonant frequencies, weakly coupled to each other 
(say, by using weak strings between neighboring pendulums), with each pendulum driven in- 
dependently by an external feed (say, via EM solenoids) with the same phase for all of them. 
By tuning the frequency of the driving feed around the resonant frequency of the pendulums, 
one may expect to observe stratified excitation of the pendulums, from long-wave, ferromag- 
neticlike strata, to short-wave, antiferromagneticlike ones. One can also expect to observe 
"magic" numbers and related effect of non-moving pendulums surrounded by strongly ex- 
cited ones. (It is worth noting that we are talking here about a cw motion, in contrast to 
the well known effect in two coupled non-driven pendulums, whereby the pendulums peri- 
odically alternate periods of zero-excitation in one and strong excitation in the other one. 
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with the frequency of the alternation being inversely proportional to the coupling.) The 
major effect here would be again due to locsitons, and the major critical condition for their 
excitation would be similar to the condition on the strength of interaction, Qa, (12.101) . to 
exceed a critical one. 

The driven pendulum array would be a great classroom-demonstration tool. However, 
although purely mechanical, it might take some effort to implement, due to the need of 
an independent feed to each pendulum. From that point of view, a much simpler, easily 
implementable, controllable, recordable, and versatile system could be an electronic array 
of individual resonant circuits, electronically coupled to each other. A long while ago, 
similar systems have been used as transmission lines, whereas a finite set of circuits has 
also been of interest for such applications as bandpass filters. However, the detailed picture 
of behavior of individual circuits inside such systems have apparently not attracted too 
much attention, with a few reasons for that. The main difference between these systems 
and the ones proposed by us, is that we need an independent feed with the same phase 
for each individual circuit, which may be arranged, by e.g. using individual cable for each 
one of them. The coupling between the individual circuits can be engineered in such a way 
that one can arrange strictly near-neighbor interaction, two-near-neighbors interaction, etc. 
Besides, the contribution of each neighboring circuit can be independently controlled, so the 
term can now be changed to any desirable function of the neighbor spacing. 

A close example, which might have important implications for large radio-frequency an- 
tenna arrays, for e. g. radio-astronomy applications, as well as in multi-dish radar systems, 
is the interaction of radiators in those arrays, if the strength of this interaction exceeds some 
critical value. 

Since in the case of electronic circuits the polarization of the feed and the spatial 
anisotropy as in (12.21) is not a factor anymore, the dipole-like interaction is simplified, as 
its spatial polarization form-factor (as in Section IIIip can be now dropped, and the 
equations of motion (I3.ip - (I3.4I) can be simplified. Furthermore, since now the direction of a 
"dipole" with respect to the "incident polarization" is not a factor, one can arrange a loop or 
ring of those circuits, instead of a linear ID array, which allows for periodic boundary con- 
dition instead of zero boundary conditions as e.g. in (13. ip . This would greatly simplify the 
theory and comparison with the experiment, on the one hand, and allow for new interesting 
effects in ring arrays (such as e. g. greatly enhanced Sagnac effect and related gyroscope 
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applications), to be discussed by us elsewhere, on the other hand. 

An interesting and exotic opportunity with ID arrays, and especially ring arrays, is a 
possibility to develop a toy model of "discrete space quantum mechanics" , whereby a wave 
function is replaced by a set of oscillating fields at discrete spatial points, instead of a regular 
wave distributed in space. This discrete-space quantum mechanics (QM) could be of interest 
for the theory of certain systems, as well as from fundamental point of view; for example, 
the ring arrays then would allow to build a theory of a Bohr-like "discrete-QM H atom", 
with finite number of primary quantum numbers, unlike an infinite numbers of quantum 
levels as in a regular-QM H atom. 

A very interesting analogy, especially from the application point of view, is magnetic 
spin resonance, in particular nuclear magnetic resonance (NMR) 19j, as well as magnetic 
resonances in finite spin systems, in particular so called molecular magnets They have 
some common points with the interacting systems considered here, in particular, two-level 
nature of resonances in both systems. 

In optical domain, an observation of the effect discussed here can be done in a few ways. 
Nanostrata and locsitons can be observed either via size-related resonances in scattering of 



laser radiation, or via X-ray or electron energy loss spectroscopy o 



The effect has promising potential for molecular computers 2l| and nanodevices. The 



the strata. 



major advantage of locsitons vs electrons in semiconductors is that they are not based on 
electric current or charge transfer. This may allow for a drastic reduction of the size limit 
for computer logic elements currently based on metal-oxide-semiconductor technology, which 
may suffer from many irreparable problems on a scale below 10 nm. As such, locsiton-based 
devices could be an interesting entry into the field, as complimentary or alternative to 
emerging technologies like plasmonics [22] or spintronics [23^]. They can offer both passive 
(e.g. transmission lines and delays), and active elements, e.g. for switching and logics. A 
ring-array may be used as a basis for a Sagnac-locsiton-based gyroscope; low locsiton velocity 
may allow for a high sensitivity in a small ring. 

Another promising application of locsitons could be biosensing devices, where target- 
specific receptor molecules either form a locsiton-supporting lattice or are attached to its 
sites; a localized locsiton occurs whenever a target biomolecule attaches to a receptor. 

Finally, exciting opportunities exist in atomic arrays and lattices with inverse population 
created by an appropriate (e. g. optical) pumping, which may lead to a laser-like locsiton 
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stimulated emitter ("locster"), to be discussed elsewhere. 



XI. CONCLUSIONS 

In conclusion, in this study of strong stratification of local field and dipole polarization 
in finite groups of atoms predicated by us earher we developed a detailed theory of the 
phenomenon in one-dimensional arrays of atom or resonant particles. In strong departure 
from Lorentz-Lorenz theory, the spatial period of those strata may become much shorter 
than the incident wavelength. By exploring nanoscale elementary excitations, locsitons, 
and resulting size-related resonances and large field enhancement in finite arrays of atoms, 
we showed that their spatial spectrum has both long waves, reminiscent of ferromagnetic 
domains, and super-short waves corresponding to the counter-oscillating neighboring polar- 
izations, reminiscent of antiferromagnetic spins. The system also exhibits "hybrid" modes 
of excitation that have no counterpart in magnetic ordering, and are more representative 
of the effect. Our theory that goes beyond Ising-like near-neighbor approximation and de- 
scribes the excitation whereby each atom interacts with all the other atoms in the array, 
reveals the existence of infinite spectrum of "exponential", or "evanescent" eigenmodes in 
the such arrays. We explored the phenomenon of "magic" numbers of atoms in an array, 
whereby resonant local-field suppression can be canceled for certain atoms in an array. We 
also demonstrated the existence of nonlinearly induced optical bistability and hysteresis in 
the system. We discussed a stratification effect similar to that in atomic arrays, which 
may exist in broad variety of self-interacting systems, from mechanical (pendulums) to elec- 
tronic circuits, to radar arrays, and to the nuclear magnetic resonance. We pointed out 
a few potential applications of the atomic and similar arrays in such diverse field as low- 
losses nano-elements for optical computers, small-size Sagnac-effect-based gyroscopes, and 
bio-censors. 
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APPENDIX A: MATHEMATICAL ASPECTS OF THE DISPERSION EQUA- 
TION dUI]) 



In this Appendix we consider mathematical properties of the Fourier series in the disper- 
sion relation (14. ip . i.e. the function 

n=l 

which is a generalized form of (14.11) [in (14.11) . p = 3 due to the chosen geometry of the 

problem - ID array of point dipoles], but we will restrict ourselves to natural numbers p. 

Note, for example, that in the case of infinite "dipole strings", parallel to each other and 

equidistantly arranged in a 2D plane, p = 2, whereas for "dipole planes" parallel to each 

other and equidistantly arranged in the 3D space, p = 1. 

Our main purpose here is to find a closed finite analytical form of the function Sp(g) 

that originated the Fourier series ( lAip . and when it is impossible in terms of more or 

less regular analytical functions, at least find a closed finite derivative dJ^T^ p{q) / dq^ with 

minimum derivative order m, which will help us to analyze in great detail the behavior of the 

function Sp(g) in its physically interesting points, e. g. at g = 0, vr, the ratio of its min/max 

values Sp(7r)/Sp(0), etc. Those derivatives are actually the main tool in our search. Indeed, 

differentiating Sp(g) in ( lAll) p — 1 times, one obtains 

d^'^^piq) _ , . p/2 sin(ng) 

dq^-^ ~ ^"^^ ^^^^ 

n=l 

if p = 2,4,6,... 

and 



(A3) 



dqP ^ n 

= (-l)(p-i)/2s^(A;), if p = 1,3,5,... 

Recalling that sin^ = (e*^ — c.c.)/2i and cos^ = (e*^ + c.c.)/2, we recognize the sums 
Yl^=i e-'^''^'' /n in (IA2p and flA3P as Taylor expansions of — ln(l — e^'"'). 

In the case of even numbers p, the sum ( IA2I) in the interval 0, 27r (which is of the most 
interest) yields: 

{7rsign((7) — o ^ , , , 

2 ' l^|e(0^2vr), ^^^^ 

0, if |g| = 0,27r. 
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It is discontinuous at g = 0, ±27r. Even numbers p are thus a "lucky" case; the result (lA4p 
is easily integrated to restore function In a physically interesting case, p = 2, by 

integrating (lA4p once we have 

^^ = 5(l^M!^, (A5) 
^ — ^ 12 

n=l 

if g G [-27r,27r], 

which is a continuous (but not smooth) function with £2(0) = 7r^/6 and S2(7r) = — 7r^/12, 

so that (S2)mm/(S2)max = -1/2. 

As to the closed integrability, one is not as lucky with odd numbers p; the summation of 
[3]) yields: 

^ ^ cosM ^ _ sin(g/2)|], (A6) 

n=l 

i.e. — s> cxD at q = 217111, 



which cannot be integrated in simple known analytical functions, but it gives a good ana- 
lytical tool to analyze the behavior of Sp(g). In the case of most interest to us, p = 3, we 
have that near the maximum of S3 at small wavenumbers 

S3(g)^^ + y (ln|g|-0 at |g| « 1, (A7) 

where S = YlTLiJ^^ ~ 1.202057, whereas near the minimum of S3, i. e. near g = vr, we have 



3S {q 



TT 



i2 



S3(g)~-— + ^^^Y^ln2, if |g-7r|«l. (A8) 

Notice that here {T,^)niin/ {^z)raa.^ = —3/4, and (fT^^/dq^ = at g = 2rmx ± tt/S. 

For the magic numbers in the case of full interaction of individual atoms with all the 
other atoms in the array, it is important to know zeros of the function S3(g). As it has been 
mentioned before [see (17. 2p and related text], it turns out that the value of the ratio g'/vr 
for the first positive root of the equation S3(g) = is very close to a small rational number: 

- = (1 + A') ■ — , with A' ^ 2.6 X 10~^ (A9) 

TT 13 

For all practical purposes, a good approximation for S3(g) is provided by: 

Sf*^(g) = cosg 



, ln[(g + Ag)|sin(g/2)|] sin^(g/2) ) 

+ HSTA5) 4 /' ^^^^^ 
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where AS = 0.01472. It coincides with the results provided by numerical summation of 
(lAll) with p = 3 with the precision better than 0.6% of S, and their zeros coincide with each 
other and with (lA9p with precision better than 10"^. 

Let us also prove the relation for (Sp)inin/(Sp)max used in f l4.4p . Indeed, 

oo 

D{tI,P) _ (Sp)min _ n=l 



n=l 



n=l m=l 



E 

n=l 

oo 



n=l 

oo 



n=l 

Note that (lAlip is valid for any number p > 1. 



n=l n=l -| . 1 



(All) 
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